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Minimal energy packings and collapse of sticky tangent hard-sphere polymers 
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We enumerate all minimal energy packings (MEPs) for small single linear and ring polymers 
composed of spherical monomers with contact attractions and hard-core repulsions, and compare 
them to corresponding results for monomer packings. We define and identify "dividing surfaces" in 
polymer packings, which reduce the number of arrangements that satisfy hard-sphere and covalent 
bond constraints. Compared to monomer MEPs, polymer MEPs favor intermediate structural 
symmetry over high and low symmetries. We also examine the packing-preparation dependence for 
longer single chains using molecular dynamics simulations. For slow temperature quenches, chains 
form crystallites with close-packed cores. As quench rate increases, the core size decreases and 
the exterior becomes more disordered. By examining the contact number, we connect suppression 
of crystallization to the onset of isostaticity in disordered packings. These studies represent a 
significant step forward in our ability to predict how the structural and mechanical properties of 
compact polymers depend on collapse dynamics. 

PACS numbers: 61.46.Bc,64.70.km,82.70.Dd,64.60.Cn 
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Over the past several decades significant research ac- 
tivity has focused on understanding dense packings of 
hard spheres, since they serve as model systems for 
atomic and colloidal liquids and glasses, jammed granu- 
lar media, and compressed foams and emulsions. An in- 
triguing property of hard-sphere systems is that they can 
be prepared in crystalline, partially ordered, and amor- 
phous packings [1[. Packings of 'sticky' hard spheres with 
contact attractions have been used to investigate self- 
assembly of colloidal particles with depletion attractions. 
Arkus et al. recently combined graph theory and geo- 
metrical techniques [2|, |3| to enumerate minimal energy 
packings (MEPs), i.e. those with the maximum number 
of contacts, for iV < 10 sticky hard spheres. Their predic- 
tions agreed with experiments on attractive colloids [J] . 

However, there have been few studies of packings of 
sticky tangent hard-sphere polymers, which can model 
polymer collapse, protein folding, and protein interac- 
tions [5|. Recent simulations [y, Q and experiments 
[8| have investigated polymer packings; however, they 
considered non-sticky spheres with only hard-core repul- 
sions, where free volume, not energy, is relevant. Thus, 
there is little understanding of how covalent bond and 
chain uncrossability constraints affect structural and me- 
chanical properties of sticky hard-sphere polymer pack- 
ings and the probabilities with which these occur. 

In this Letter, we perform exact enumeration studies 
of MEPs for sticky, tangent, monodisperse hard-sphere 
polymers (both linear and cyclic) and contrast the re- 
sults with those for sticky hard spheres without polymer 
constraints. Our studies begin to address several overar- 
ching questions: 1) How do the probabilities for obtaining 
polymer MEPs differ from those for sticky hard-sphere 
MEPs? and 2) How do the properties of single compact 
polymers depend on collapse dynamics, e.g. do they col- 
lapse into crystalline or amorphous clusters? 



Our results show that polymer constraints reduce 
the ways in which hard spheres can be arranged into 
MEPs, and the strength of this effect varies for differ- 
ent macrostates (i.e. structurally distinct packings). We 
demonstrate that the large reduction in the number of 
arrangements may be understood in terms of dividing 
surfaces. These split polymer packings into disjoint re- 
gions and eliminate particle-label permutations that do 
not correspond to polymer chains. We find that polymer 
MEPs with intermediate structural symmetry are more 
frequent relative to the monomer case, where entropy fa- 
vors low symmetry MEPs [4|. 

In addition, using molecular dynamics (MD) simula- 
tions of temperature quenches at various rates T, we show 
that single chains display glassy dynamics during col- 
lapse, and that the final polymer packings depend on T. 
In the slow quench rate limit, the chains undergo a sharp 
[Oj transition to crystallites, with a jump in the energy 
and number of contacts N^ (including covalent bonds) at 
temperature T = Tmeit- The crystallites possess a close- 
packed core surrounded by a "surface" whose size and 
disorder increase with \T\. For slow quenches, A''c at Tmeit 
jumps from below the minimal number A^™'" = SN — 6 
required for mechanical stability [lOj to Nf""", where a 
significant fraction of the monomers possess 12 contacts. 
In the large \T\ limit, the clusters are disordered with 
< ^mm contacts even as T ^' 0, showing that rigidifica- 
tion can hinder crystallization. 

We first describe exact enumeration methods for 
monomer and polymer MEPs [20|. To generate possi- 
ble packings for a given number of spheres N and contact 
number Nc, we iterate over all A^ x A^ adjacency matrices 
A satisfying J2i>i A-ij = Nc. The elements of A are 1 for 
contacting particles, and for non-contacting particles 
and diagonal entries. Covalent bonds link sticky spheres 
to form a polymer chain with length A'^; Ai^i^i = 1 for 



TABLE I: Statistics for MEPs with N spheres and Nc contacts. M is the number of macrostates, fr, fp, and fm are the 
fraction of microstates obeying minimal rigidity constraints that also satisfy hard-sphere constraints, respectively for rings, 
linear polymers, and monomers, and f2r, Qp, and Qm are the total numbers of microstates satisfying both minimal rigidity and 
hard-sphere constraints. Values for f and Q do not account for chiral twins [2|. In agreement with [2], we find 1 and 4 floppy 
macrostates (in the fc — ^ oo limitlM), respectively for N = 9 and (A'^ = 10, Nc = 24). However, we find 2 and 55 more rigid 
macrostates*'* for these cases [13, 1201 . Adjacency matrices and coordinate solutions for all microstates are available online 
\2M. — indicates data not available. 
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1 < z < A^ for linear chains, and additionally Ai^n = 1 
for rings. The distinction between permanent covalcnt 
and thermally fluctuating noncovalent bonds is not im- 
portant for static packings; we include both types in Nc- 
We enumerate all adjacency matrices satisfying the 
above conditions and then identify those that also ful- 
fill hard-sphere and minimal rigidity constraints. Hard- 
sphere constraints imply that the center-to-center dis- 
tances r,j between unit spheres i and j obey r^ > 1, 
where the equality holds for contacting pairs. Necessary 
conditions for rigidity are that each monomer possesses 
at least three contacts and Nc 



13 



> A^f " 

To enforce these constraints, we irnplemented geomet- 
rical rules developed by Arkus et al. [3, |3| that eliminate 
invalid adjacency matrices. For the remaining configura- 
tions, we solved the system of quadratic equations 



= di 



(1) 



for sphere positions Vi to an accuracy of 10~^. We also 
calculated the dynamical matrix (all second derivatives of 
the energy in Eq. ^ with respect to monomer displace- 
ments) for all configurations, which allowed us to identify 
rigid (with 3A^ — 6 nonzero eigenvalues) and floppy con- 
figurations |13| (with fewer nonzero eigenvalues) [20|. 

From this procedure, we obtain microstates and 
macrostates for a given A^ and Nc that satisfy hard- 
sphere and minimal rigidity constraints and the relevant 
polymeric constraints. Each macrostate is characterized 
by an adjacency matrix that is nonisomorphic to and a 
set of interparticle distances {rij^ that is different from 
those characterizing other macrostates [20[. With this 
definition, no macrostate can be rotated or reflected such 
that it yields a different macrostate. Every connected 
sticky hard-sphere macrostate admits a linear polymer 
macrostate [14| . Thus, sticky-sphere and linear polymer 
packings have identical macrostates. We have also veri- 
fied this for ring packings for A^ < 10. 



A microstate is a particular labeling of the parti- 
cles 1 through A^ that comprise a A^-particle macrostate 
with A'c contacts. Many microstates correspond to each 
macrostate due to particle permutations for monomer 
packings [j|, and for polymers, the multiple possible 
paths through a given macrostate. The total number 
of microstates flm, ^p, and Q,r is given by the sum 
of microstates for each macrostate for monomers, lin- 
ear polymers, and rings, respectively [20|. For monomer 
packings, which lack covalent bonds, the number of mi- 
crostates for each macrostate (ignoring chirality) is given 
simply by a geometric factor fi^ = P^, where Pi is the 
number of allowed permutations of particle indices for 
macrostate i [3|. For polymer packings, the number of 
microstates is not given by this relation since one must 
ensure that particle indices are consecutive. 

Exact enumeration results are displayed in Table U 
which shows the number of macrostates M, fraction / 
of adjacency matrices with Nc contacts obeying minimal 
rigidity that also satisfy hard-sphere constraints, and r2,„, 
57p, and O^ for 5 < A^ < 10. /is the probability to obtain 
a packing for an 'ideal' protocol that samples adjacency 
matrices uniformly. From TableU we see that / decreases 
approximately exponentially with A^ for N > 5, and even 
faster for A^ > 9. Part of the reason for the strong de- 
crease in / between A^ = 9 and 10 is the decrease in 
macrostates from 52 to 3. This occurs because A^ = 10 
MEPs possess Nc = A"™'" 4- 1, which exceeds the number 
of degrees of freedom. Eq. [T]is then overconstrained, and 
its solutions possess special symmetries. The increase in 
Nc signals the onset of crystal nucleation, and the for- 
mation of a close-packed core. The ability to enumerate 
the numbers of isostatic {Nc = A"""") and hyperstatic 
{Nc > A"""") packings will yield insight into systems 
where glass and crystallization transitions compete. 

For the A^ studied here, hard-sphere constraints are 
more difficult to satisfy for minimally rigid polymer pack- 



ings compared to monomer packings: fr<fp< fm 115 
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FIG. 1: (a) Schematic of dividing surfaces S and S' (colored 
triangles formed by monomers {B, C, D) and [D, E, F), re- 
spectively) in a macrostate for A^ = 8. For S, region J consists 
of monomer A and region K of monomers {E, F, G, H), or 
vice versa, (b) Fraction of microstates for packings from each 
symmetry group for cyclic (open circles) and linear (down- 
ward triangles) polymers, and monomers (filled circles) with 
N = 8. Results in (b) do not account for chiral structures. 

A key mechanism for the reduction in / is the occur- 
rence of "dividing surfaces" in polymer packings. A di- 
viding surface is any minimal subset of a connected clus- 
ter of contacting monomers that geometrically splits it 
into two. Any polymer path that traverses a dividing 
surface that does not also topologically divide the poly- 
mer is blocked and invalid. Specifically, if m consecutive 
monomers i + 1^ . . . ,i+m occupy an ?7i- monomer dividing 
surface S, any polymer path where the sets of monomers 
J and K divided by 5* are anything other than 1, 2, . . . , i 
and i + m + l,i + m + 2,...,N (or vice versa) is blocked. 
In other words, any path that starts in J, enters S, and 
traverses it (passes through all monomers in S) is blocked 
unless it traverses all monomers in J before entering S. 
Fig. [ija) schematically depicts the sets J and K and two 
dividing surfaces for a A^ = 8 macrostate. By definition, 
blocking does not occur in monomer packings. 

In Table IH we see that the blocking effect increases 
sharply with A^ since fr/fm and fp/.fm decrease signif- 
icantly. Blocking also reduces [15[ the fractions of al- 
lowed ring microstates relative to those for linear poly- 
mers fr/ fp since rings do not possess chain ends. An- 
other clear feature in Fig. [Ub) is that blocking changes 
the relative frequencies with which macrostatcs of differ- 



ent symmetries are populated. Ring and linear polymer 
packings are more likely to possess intermediate sym- 
metry than monomer packings, whereas the opposite is 
true for macrostates with the lowest and highest symme- 
tries. Highly symmetric macrostates possess many dis- 
tinct blocking surfaces, and low symmetry macrostates 
possess a surplus of closed trimers as shown in Fig.[lja). 

The enumeration studies illustrate an interesting com- 
petition between energy and entropy in large systems. 
For N > 10, MEPs are overconstrained with Nc > N^'"". 
This suggests that if the system becomes trapped in a 
metastable state {e.g. with Nc — A"™'"), rearrangements 
into MEPs will be slow because of their low entropy. 
Thus, glassy dynamics in single polymer chains should 
be observable in systems quenched at varying rates. For 
ksT :§> |e|, where — e is the contact energy, polymers will 
adopt random-coil configurations with Nc ^ A™'". As 
the polymer is cooled, one expects quench rate effects to 
become important when Ac ~ A"™'" [16[. 

To demonstrate glassy dynamics for single linear 
polymer chains, we employ MD simulations in which 
monomers interact via the potential energy 
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where k is the spring constant and Z) = 1 is the 
monomer diameter. The temperature T is controlled 
via a Langevin thermostat. The unit of time is t = 
■y/m-D^/e, where m is the monomer mass. The cutoff 
radius rd D = cxd for covalently bonded monomers and 
1 -|- y2e/fc for noncovalently bonded monomers. C/harm 
reduces to the energy for tangent sticky hard spheres 17 1 
in the limit fc — > cxd and possesses the same MEPs. For 
A^ < 10, the MEPs from simulations agree with those 
from complete enumeration for k > 1600e [tc ^ 1.04Z)). 

Figure [5Ja) shows the potential energy per particle 
—U/Ne for different quench rates T. At low T, a sharp 
transition between coils and crystallites [9|,|ll| is observed 
at Tnioit — O.STe/fcs. The crystallites consist of a close- 
packed core with Ncp monomers (each with 12 contacts) 
and a less-ordered exterior. The crystallization transi- 
tion corresponds (Fig.[2Ib)) to a sharp transition in A"cp, 
which implies a change of symmetry within the core, from 
liquid-like to close-packed. 

At higher rates, the dynamics becomes glassy near 
Tmoitj and the systems do not approach the ground state 
energy even as T — )- 0. We associate the suppression of 
crystallization with the onset of rigidity. Evidence for 
this is given in Fig. [2][c). The data show two "critical" 
quench rates: t* and t** . For \f\ < \f*\, the jump 
in Ac and Ncp resembles a first-order transition. For 
|r| > |r**|, the systems do not form minimally rigid 
clusters even at T = 0. Even though the critical rates 
and Tmcit are A^-dependent, the trends are clear. For 
A ~ 100 systems, we estimate \kBT* /e\ ~ 10^'^/r and 
\kBt**/e\ ^ 10-3/t. 
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FIG. 2: (a) Potential energy per particle { — U/Ne) (b) num- 
ber of particles with 12 contacts (Ncp), and (c) total num- 
ber of contacts (Nc) versus ksT/e for single linear polymers 
with A'^ — 100 at different quench rates. Data (top to bot- 
tom, panels b-c; bottom to top, panel a) are for quench rates 
ksTT/e = -10"^ -10"*, -10"^ -10"'^ and -10"''. The 
critical quench rates are IksT* /e\ ~ 10"'/t and |fcsT**/e| ~ 
10"^/r. All results are averaged over several independent 
initial configurations. The horizontal (vertical) dotted lines 
indicate Nc = A^""" {kBT/e = 0.37). 

The effects of quench rate on end states of quenches 
to T = are visualized in Fig. [31 Monomers are color- 
coded by the number of contacts; dark blue (red) in- 
dicates close packing (^ 12 contacts). The left panel 
shows a typical configuration after a fast quench with 
ksT = — 10~'*e/r; wc sec a small close-packed core sur- 
rounded by a disordered exterior. The middle and right 
panels show a collapsed structure at T = from a slow 
quench {ksT = — 10"^e/r). The close-packed core is 
much larger, and the exterior is more crystalline. The 
large gaps visible in the rightmost panel indicate the or- 
der is hep, and the structure is stack- faulted p^ . 

We examined minimal energy packings of sticky tan- 
gent hard-sphere linear and cyclic polymers, and com- 



FIG. 3: Collapsed structures at T = for a single N = 100 
linear polymer using two quench rates: kBTr/e = — 10~* 
(left) and -10"'' (middle, right). The packing in the right 
panel is rotated compared to that in the middle panel to show 
its hexagonal planes. 

pared them to monomer packings for small N. The 
packings are the same, but polymer packings possess sig- 
nificantly smaller entropies compared to monomer pack- 
ings due to dividing surfaces, which arise from covalent- 
bond constraints. Entropic suppression via blocking is 
strongest for structures of both very high and low sym- 
metry. In both monomer and polymer cases, the fraction 
of states satisfying hard-sphere constraints decreases at 
least exponentially with increasing N, and faster when 
Nc > iV""". We also performed MD simulations of single 
linear chains with larger N, which link glassy dynamics 
to the onset of rigidity. This work sets the stage for fu- 
ture studies that investigate whether cooperative dynam- 
ics from chain connectivity and uncrossability constraints 
improves or impedes glass-forming ability of single poly- 
mers compared to colloidal systems. 

We thank V. N. Manoharan for helpful discussions. 
Our results were obtained using the Boost Graph Library, 
a modified version of N. Arkus' structure solver [3|, and 
LAMMPS f[^. Support from NSF Award No. DMR- 
0835742 and an Anderson Fellowship from Yale Univer- 
sity is gratefully acknowledged. 



SUPPLEMENTARY MATERIAL FOR "MINIMAL 

ENERGY PACKINGS AND COLLAPSE OF 

STICKY TANGENT HARD-SPHERE 

POLYMERS" 

In this supplementary material, we provide additional 
details concerning the methods employed to: (1) enumer- 
ate exactly all micro- and macrostates and (2) assess the 
rigidity for monomer and polymer packings composed of 
sticky, monodisperse tangent hard spheres. 



Exact enumeration method 

The exact enumeration method consists of several steps 
including looping over all adjacency matrices satisfying 
the appropriate constraints for monomer and polymer 
packings, identifying those adjacency matrices that sat- 
isfy hard-sphere constraints, and then solving for their 
Euclidean positions. 

The adjacency matrix A for a iV-particle system is a 



N X N symmetric matrix whose elements are 1 for con- 
tacting particles and for noncontacting particles. For 
monodisperse hard spheres, 



A, = 1 if 



= D, 



A, 



if ry > D, 



(3) 



where r^ = |ri — fj 



is the center-to-center distance be- 
tween spheres i and j and D is their diameter. By con- 
vention, the diagonal entries satisfy An = 0. 

The number of permutations for symmetric N x N 
matrices with N const constrained elements {e.g. covalent 
bonds) is 



Pc. 



[{N'^ -N)/2]\ 



Afco„«t![(A^2-iV)/2-iVco„«t]!' 



(4) 



where the {N'^~N)/2 terms arise because we only need to 
consider entries above the diagonal since A is symmetric. 
The number of permutations Pcont for A with N^ sticky- 
sphere contacts has the same form as Eq. |4l 



[{N'^ -N)/2]\ 



N,}.{{N^-N)/2-Nc)V 



(5) 



For the packings considered in this study, the con- 
strained elements correspond to covalent bonds, which 
are fixed to be 1 not 0. Thus, the number of permuta- 
tions is 



^cct(-^const) 



-Wconst 



(Wc-Wco„.t)! [^^^-(Wc-Wconst)] ! ' 



(6) 



where A'const ~ 0, iV — 1, and N , for monomers, linear 
polymers, and rings, respectively. Specifically, for linear 
and ring polymers, Ai^+i = 1 for \ < i < N . Rings 
must also satisfy the constraint Ai^jsi — \. In contrast, 
monomer packings do not have explicitly constrained off- 
diagonal elements. 

Thus, the numbers of adjacency matrices for 
monomers, linear polymers, and rings [23 for TV particles 
and Nc contacts (including covalent bonds for polymers) 

are Irn. — -front , 



Pir>= -Pcct(iV-l), 



r JV^-3JV+2 i| 



{Nc -{N- l))!(iv!^M+2 _ (TV, -{N- 1))! ' 



and 



P, =Fcct(iV), 



[(iV!_3^], 



(7) 



(8) 



{Nc-N)\{^^^^-{Nc-N))\ 



We then loop through all permutations Pm, Pip, and 
Pr of A for a given N and Nc. Since the entries of A 



are ones and zeros, different adjacency matrices corre- 
spond to unique binary numbers {i.e. each microstate 
corresponds to a particular adjacency matrix and unique 
binary number). We enumerate all binary numbers using 
sequential binary permutations (from the C-l — h Standard 
Template Library nextjpermutation() function). Note 
that considering polymers leads to an exponential reduc- 
tion in the effort required for exact enumeration. 

Since we are interested in minimal energy packings 
(MEPs) — those with the maximum number of contacts 
Nc, we focus on packings with Nc > iV™'", where 
^mm ^ 2^ — 6 is the minimal number of contacts re- 
quired for rigidity. For 4 < A^ < 9, we verified that no 
microstates with Nc > TV""" satisfy hard sphere con- 
straints (as shown previously [2|), and for A^ = 10 no 
microstates with Nc = N"^™ + 2 exist, i.e. MEPs for 
4 < TV < 9 possess TV^ = TV™'" and for TV = 10 possess 
TVc = TV™ + l. 

To eliminate adjacency matrices that do not satisfy 
hard-sphere constraints, we employed the complete set 
of geometrical rules for the adjacency matrices for TV < 7 
provided in Ref. [3|. Specifically, we implemented rules 
1-12, 14—18 outlined on pages 295-318. However, we did 
not employ the triangular bipyramid rule (discussed on 
pages 40-48 of [Sj]), neither the version for iterative pack- 
ings nor that applied to new seeds. Instead, for TV > 7 we 
inserted those minimally rigid adjacency matrices not re- 
jected by the geometrical rules into a modified version of 
Arkus' Euclidean structure solver. The structure solver 
makes a random initial guess for particle coordinates and 
then uses Newton's method to solve the contact equa- 
tions implied by the adjacency matrix, while enforcing 
hard sphere constraints. We checked for convergence of 
the structure solver by increasing the maximum number 
of 'initial guesses' for the coordinates; 3TV'^ initial guesses 
arc sufficient to solve all structures to an accuracy in po- 
sitions of 10~^. 

Closely associated with A is the distance matrix D 
whose elements are Dij = r.ij. For TV < 10 and TV = 10, 
TVc = 25, nonisomorphic A (identified using the Boost 
Graph Library's isomorphism() function) correspond to 
different macrostates because the elements of D are dif- 
ferent ^. However, for TV = 10, TV^ = 24 graph 
(non) isomorphism (as given by isomorphism ()) is insuffi- 
cient to completely distinguish macrostates. While there 
arc 286 nonisomorphic graphs, 8 of these produce coor- 
dinate solutions that are identical to those produced by 
other graphs. All eight of these correspond to 'switching' 
a noncovalent bond in such a way that the same coordi- 
nate solution is produced. This reduces the total number 
of macrostates to 278. 

For example, graphs 128 and 158 produce the same 



set of coordinates [2l| , but differ in that the former has 



a noncovalent bond between particles 1 and 6, while the 
latter possesses a noncovalent bond between particles 4 
and 7. We assume that graphs 128 and 158 correspond to 



the same macrostate because they possess the same D, 
but (since the particles are distinguishable) the coordi- 
nate solutions correspond to different microstates |24| . 
We therefore have retained aU microstates for all 286 
graphs, but assign them to 278 macrostates. Adjacency 
matrices for graphs 128 and 158 and their common coor- 
dinate solution are shown below in ([9]); red entries indi- 
cate the 'switched' bond. 



1110 11111 
1010110001 
1101100011 
1010110100 
1110 10 
1101101000 
10 10 10 
10 10 10 10 
10 10 10 1 
1110 10 



Graph 128: 



Graph 158: 



1110 1111 
1010110001 
110 110 11 
1010111100 
1110 10 
10 110 10 
10 10 10 10 
10 10 10 10 
10 10 10 1 
1110 10 
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following harmonic intcrparticle potential 






(10) 



where D = 1 and k must be large enough so that no 
'2nd-nearcst neighbors' interact. We have examined the 
number of rigid macrostates as a function of k. For suffi- 
ciently high precision in the coordinates (one part in 10^ 
or better), all macrostates for 4 < TV < 8 are rigid {i.e. 
the dynamical matrix possesses 3iV — 6 positive eigen- 
values) for fc/e > 10. For the eigenvalue threshold, we 
assumed that eigenvalues e^ > 10~^ were nonzero and 
positive. For iV = 9, we find one 'floppy' macrostate for 
all fc/e > 400, consistent with Ref. [2[. We also verified 
that all N = 10, Nc = 25 macrostates possess are rigid. 

For TV = 10, Nc = 24, however, the results show non- 
trivial dependence on the numerical precision of the co- 
ordinate solutions. When we solve for the coordinates 
to a precision of one part in 10^, we find 4 nonrigid 
(floppy) macrostates (in agreement with Ref. [2[) for all 
k/e > 400. However, when the coordinates are solved 
to a precision of only one part in 10®, one of the floppy 
configurations becomes 'rigid' in the same range of k due 
to insufficient precision. 

Figure m shows the eigenvalue spectra for the A^ = 10, 
Nc = 24 macrostates for a precision of one part in 10^ in 
the coordinate solutions. The eigenvalues are displayed 
from smallest to largest: j = I to 37V. The 3iV — 6 pos- 
itive eigenvalues for rigid macrostates are well-separated 
from the 6 eigenvalues that correspond to rigid transla- 
tions and rotations as shown in panel (a). The eigen- 
values corresponding to rigid translations and rotations 
(indexes 1-6) are zero to within our numerical precision. 
For floppy macrostates, the floppy eigenvalues are also or- 
ders of magnitude below those that correspond to finite- 
energy normal modes. 

Since nonrigid macrostates have fioppy modes that 
can be activated with no energy cost, each nonrigid 
macrostate possesses a continuum of coordinate solu- 
tions. However, we have verified that none of the four 
floppy macrostates for iV = 10, A^c = 24 can be trans- 
formed continuously into one another (without increas- 
ing C/harm)- Specifically, we have shown that different 
coordinate solutions corresponding to the same floppy 
macrostate differ by ^ 10~'^D, while transitions between 
floppy macrostates would require displacements ~ D. 



(9) 



Rigidity analysis 

The rigidity of all macrostates generated by our exact 
enumeration algorithm was assessed by calculating the 
eigenvalues of the dynamical matrix [25|, assuming the 
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sets of components 
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for A'' < 10, but we anticipate it will become increasingly 

important for enumeration of MEPs at higher A'^. 

To see this, note that if the particles are distinguishable, 

we can assign 'colors' to them. Suppose 1 = "red", 4 = 

"green", 6 — "blue", and 7 = "yellow". Thus 128 contains 

a red-blue bond while 158 contains a green-yellow bond, 

and these are distinguishable at the microstate level. 
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